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ABSTRACT 



This work addresses the design and analysis of the discrete longitudinal 
autopilot for application to the bank-to-turn (BTT) missiles. 

In the development the classical design and analysis were reviewed 
using the continuous uncoupled pitch channel autopilot for circular airframe. 
Then,app1ying analog-to-digital conversion, the corresponding discrete 
autopilot was designed. The performance of both continuous and discrete 
open loop systems was analyzed according to the desired requirements. 

In the following section, utilizing modern control design techniques for 
the discrete pitch channel autopilot, the control law was designed, 
assuming availability of all the states and taking the desired poles to be at 
the same place as in the classical design. Therefore the discrete state- 
feedback autopilot was obtained and analyzed. The next step was the design 
of an estimator which estimates the entire state vector, given 
measurements of the portion of the state. The performance of the estimator 
was analyzed too. 

Finally coupling the discrete pitch and roll channel autopilots the 
overall performance of the discrete system was obtained and analyzed. The 
resulting design was also proved to be robust. 
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TABLES OF SYMBOLS AND ABBREVIATIONS 



BTT Bank-to-T urn 

CBTT Coordinated Bank-to-Tum, minimum sideslip, positive a, 
9# i 180 deg. 

C] rolling moment coefficient 

1^ pitching moment coefficient 

Cm* slope of curve of pitching moment coefficient C m vsa 

Cmsp change in Cm per degree pitch control incidence, 6 p 

C N normal force coefficient 

Cffc slope of curve of normal force coefficient C N vsa 

Cm£ p change in Cf< per degree pitch control incidence,^ 

C Y side force coefficient 

d reference length for coefficients (= 2 ft) 

I yy moment of inertia about y^ axis 

I K moment of inertia about axis 

moment of inertia about x^ axis 
K a autopilot pitch acceleration error gain 

K Y p CBTT autopilot coordination branch gain 

p roll rate about x^ 

p roll acceleration about x^ 

POC preferred orientation control 

<f dynamic pressure 

q pitch rate about 

q pitch angular accelaration about y^ 
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constant or equilibrium pitch angular rate 
yaw angular rate about 5^ 
yaw angular acceleration about 5^ 
reference area for coefficients (=tf ft 2 ) 

Skid-to-Tum, roll attitude stabilized 

velocity component in Jfc direction 

velocity component in Ife direction, assumed to constant 

constant missile flight path velocity 

missile velocity vector 

velocity component in 5^ direction 

body-fixed roll axis, along axis of symmetry, positive forward 
body-fixed pitch axis, positive starboard 
body-fixed yaw axis, forms right handed orthogonal system 
with * and ^ 

vehicle axis in downward direction along local gravity 

vector , approximated as inertial axis 
achieved normal acceleration in Zq direction 
commanded normal acceleration in z B direction 
achieved normal acceleration in y^ direction 
achieved normal acceleration in % direction 
achieved normal acceleration in y^ direction 
normal acceleration command from guidance computer in 
F y direction plus anti-gravity bias command 
normal acceleration guidance command in % direction 
normal acceleration guidance command in y^ direction 



<pc 



<p 



<Pe 

9 

♦ 

Spc 

Sy 

S Yc 

$Rc 

a* 

a 

i 

a 

* 

* 

OPTSYS 

ORACLS 

X 



roll attitude command from guidance computer, zero degrees 

in -Zy direction and 90 degrees in direction 

roll attitude, zero degrees in -z^ direction and 90 degrees 

in £ direction 

roll attitude error, <p c -<p 

Elevation Euler Angle, second rotation, J(q cos(p-rsin(p)dt 
Azimuth Euler Angle, first rotation about y^ J(q cos<p+rsin<p)dt 
pitch control incidence (positive tail incidence produces 

negative pitching moment) 

commanded pitch control incidence, S p 

yaw control incidence (positive tail incidence produces 

negative yawing moment) 

commanded yaw control incidence, Sy 

roll control Incidence (positive tail incidence produces 

positive rolling moment) 
commanded roll control incidence , Sr 
constant or equilibrium angle-of-attack 

angle-of-attack 

angle-of-attack rate 

angle of sideslip 

sideslip angular rate 

Optimal Systems Control Program 

Optimal Regulator Algorithms for the Control of Linear 

Systems 

plant state = N s x 1 or n x 1 
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N s 

U 

No 

F 

A 

G 

8 

y 

H 

J 

K 

L 

T 



number of states 

control = N c x 1 or m x 1 

number of controls 

continuous plant system matrix 

discrete plant system matrix 

continuous plant control input matrix 

dicrete plant control input matrix 

output measurements 

plant output matrix 

plant direct transmission matrix 

control gain (law) 

estimator gain 

sampling period 
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I. INTRODUCTION 



Missile systems designed to meet the threats of the future will require 
capabilities far beyond of those currently available. Future requirements 
dictate the need for dramatic increases in Short Range Air-to-Air Missile 
(SRAAM) maneuverability and guidance accuracy. Many current missile 
designs, which are characterized by cruciform airframes and skid-to-turn 
(STT) steering, cannot be extended to meet these requirements due to low 
aerodynamic efficiency. Bank-to-turn steering provides the capability to 
design asymmetric missile airframes which optimize aerodynamic 
performance in one plane and thus obtain the SRAAM maneuverability and 
accuracy requirements. Furthermore the ramjet engine chin inlet 
configuration is believed to provide a greater range capability than other 
inlet configurations, and BBT control is required to satisfy the sideslip 
constraints imposed by chin inlets. 

In addition, to its potential advantages, BBT steering introduces some 
technical concerns which must be carefully evaluated. For example, the 
methodology for designing a bank-to-turn autopilot is not well developed. 
Such a design must take into account the aerodynamic and kinematic 
coupling terms as well as allow for operation a low signal levels (i.e, small 
angle of attack) when the preferred roll orientation is poorly defined. In 
addition, the coupling of body motion into the guidance signals is another 
major concern for BBT systems. Recent analyses have indicated that some 
skid-to-turn (STT) systems can tolerate a limited amount of radome- 
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induced instabilities without severe performance degradation [Ref. 1 ]. In a 
BBT system a coupling loop is closed through roll rate as well as pitch and 
yaw rates [Ref .2]. It is not known whether BBT systems can tolerate 
coupling induced instabilities. Another concern is the interaction of BBT 
control with missile functions such as detection (seeker), guidance signal 
processing , control surface effectiveness, etc. 

All the above concerns must be investigated before BBT steering can be 
considered a viable method to control high performance tactical missiles. 

Many missile programs were initiated during the past decade to 
improve the capability of steering a tactical missile via BTT control. The 
results have greatly advanced the understanding of the various missile 
subsystems. 

In the autopilot area, specifically, many types of autopilots have been 
found which force the missile to roll or bank so that the steering maneuver 
occurs with the missile airframe oriented in a specific or preferred 
direction with respect to the incoming airsteam. This entire class of 
autopilots referred to as Preferred Orientation Control (POC) autopilots. 

The main criterion for the selection of one particular type of autopilot 
is based upon the guidance, airframe, and propulsion system requirements. 
Generally missiles with either one or two planes of symmetry use a POC 
autopilot which forces the missile to bank in order to turn (BTT) like an 
aircraft, and if this motion is coordinated it is usually referred to as CBTT. 

In order to take full advantage of CBTT control, planar airframes have 
been designed to increase lifting capability in one direction without the 
weight and drag penalty associated with orthogonal lifting surfaces [Ref. 1 ]. 
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These airframe have aerodynamic properties which have the potential to 
enhance CBTT control. 

The present work addresses the design and analysis of the discrete 
pitch channel autopilot for application to the bank-to-turn (BTT) missiles, 
with circular airframe configuration. The circular airframe is chosen 
because it is stable for all angles of attack. 

In the development the classical design and analysis of the uncoupled 
pitch channel autopilot as were developed in reference [Ref.3], are reviewed. 
Using state-space representation a tenth order system is formulated and 
the Step Response, Bode and Pole-Zero Map plots are obtained and analyzed 
according to the desired requirements. Then the continuous system is 
converted into the corresponding discrete system using Analog to Digital 
conversion from ORACLS, and analyzed according to the Step-Response,and 
Pole-Zero Map plots obtained. A comparison of the plots of both the 
continuous and discrete systems shows no significant differences. 

Next, utilizing modem control design techniques and assuming that we 
have all the states at our disposal for feedback purposes, the control law is 
designed using the Ackermann’s formula (Ref. 6] and taking the desired poles 
to be at the same place as in the classical case. Using the control law gains 
the discrete state-feedback autopilot is designed and the Step Response,and 
Pole-Zero Map plots are obtained and analyzed. Then, since we usually know 
only a portion of the state variables, an estimator is introduced in order to 
estimate the entire state vector given measurements of that portion. The 
poles of the estimator are moved to the left of the original poles of the 
system in order to get faster response. Plots as above are obtained and an 
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analysis is made. Next the designed discrete pitch channel autopilot is 
coupled with the discrete roll channel autopilot and the whole discrete 
system is analysed and proved to be robust. Finally conclusions and 
recommendations for future study are stated. 

The analysis in all the above cases was performed using the existing at 
NPS Optimal Systems Control Program (OPTSYS) for the continuous system. 
Optimal Regulator Algorithms for the Control of Linear Systems (ORACLS) 
for the discrete system, and Pole Placement and Robustness Design Program 
(POPLAR) for the robustness of the discrete system. 



22 



II. CLASSICAL LINEAR DESIGN AND ANALYSIS OF LONGITUDINAL UNCOUPLED 
CHANNEL AUTOPILOT FOR CIRCULAR AIRFRAME 



A. GENERAL 

The Initial phase In the design of the CBTT autopilot Involved the design 
of Individual, uncoupled channels, pitch, yaw, and roll, with prescribed 
relationships between speeds of response which would meet the CBTT 
requirements when coupled. 

The design technique for the uncoupled pitch channel was classical, using 
a combination of Frequency Response and Root Locus techniques [Ref. 2], to 
achieve practical bandwidths and In turn provide the range of required 
missile body angular rates and control motions. In addition, the resulting 
design Is found to be robust. 

The whole design was done for a flight condition at M=3.95 and 
Alt1tude=60 KFT. 

The requirements for the classical design technique, as stated In 
[Ref. 2] are: 

1. Hi gh Freouencu Attenuation in Actuator Command Branch 

> 15 dB at 100 rad/sec and zero angle-of-attack. This requirement 
limits autopilot speeds of response. 

2. Relative Stabilit u: 

Gain margins > 6 dB, Phase margins * 30 degrees with a goal of 
12 dB and 50 degrees. 
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3. Acceleration Time Response 

a. 63 percent time constant of 0.5 seconds for a step command of 
acceleration at the flight condition of interest (M=3.95, A1titude=60 KFT) 
and small angle of attack. This response is representative of a tactical 
missile of this size. 

b. Overshoot < 10 percent. 

c. Zero steady state error in acceleration to reduce variations of 
guidance navigation gain. 

In this chapter the airframe configuration, the aerodynamic model and 
control law for the pitch channel are reviewed. The transfer functions for 
circular airframe which is stable for all angles of attack [Ref. 21, are 
derived. An analysis for the uncoupled pitch channel, as regards the 
acceleration response, the body angular rate and control surface deflection, 
in terms of time resposnes and pole-zero map plots for both the continuous 
and discrete system and frequency responses for the continuous system, is 
made, using the existing at NPS Optimal Systems Control Program (OPTSYS) 
for the continuous system, and Optimal Regulator Algorithms for the Control 
of Linear Systems (ORACLS) for the discrete system. Finally a comparison of 
the continuous and discrete system by analyzing the above results is made. 

A general block diagram of the BTT autopilot is shown in Figure 2.1. 

B. AIRFRAME CONFIGURATION 

An investigation has been conducted by NASA [Ref.4] to compare the 
experimental aerodynamic characteristics of a low-drag missile concept 
with a body of circular cross section to one with a body of 3:1 elliptical 
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BTT autopilot 
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Figure 2.1 Bank-to-Turn Autopilot. 



cross section, the bodies having identical cross section area distributions. 
The concepts were of mono-wing design with constant wing span. Tail 
surfaces were located flush at the body base with ±30 degrees dihedral. 
Wind-tunnel tests were performed at Mach numbers from 0.5 to 4.63 and at 
angles of attack from about -5° to +28 °. 

The comparison shows no significant subsonic normal force differences 
at low angles of attack; however, at supersonic speeds, the elliptical 
concept increasingly provides greater normal force up to Mach 2.5 to 3.0, 
beyond which an incremental increase of about 25 percent holds through the 
angle of attack range. More pronounced nonlinearities in pitching moment 
occur at subsonic speed for the elliptical concept, as well as less 
longitudinal stability at all test Mach numbers. However, levels of 
directional and lateral stability are increased, especially at the higher 
angles of attack. 

This work uses the circular airframe configuration which is shown in 
Figure 2.2 and it is taken from [Ref. 4]. The circular cross sectional body has 
a closure ratio A^/A m4X of 0.69 and the A™* occurs at 68 percent body 

length. 

C. UNCOUPLED LINEAR PITCH CHANNEL AUTOPILOT 
I. Aerodunamic Model 

The following addresses the linear design and analysis of the 
continuous uncoupled pitch channel autopilot for the circular airframe. 
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Figure 2.2 Model of Circular Configuration- 1/6 Scale. 



A linearized aerodynamic model was developed for stability studies. 
The method used is an extension of the linearization technique used for 
skid-to-turn (STT) aerodynamic models. 

The following three assumptions were made: 

a. Plane X B -Z Q of Figure 2.3 is the maneuver plane. 

b. Missile is trimmed in pitch (i.e, M y =0, at fixed values of a, q, and 

Sp). 

c. Missile roll rate is constant. 

The resulting model is shown in block diagram form in Figure 2.4. 
Aerodynamic stability derivatives for M=3.95 are provided in Table I. 
2. Pitch Control Law 

The pitch control law for the circular airframe is shown in Figure 2.5 
as given in [Ref. 21. Lag-leads were used to prevent guidance noise 
saturation problems. The rate error compensation determines the high 
frequency attenuation and was used to minimize the effect of aerodynamic 
variations on acceleration time response. The acceleration error 
compensation determines the acceleration time response. An integrator was 
used in the acceleration error branch to satisfy the guidance requirement of 
zero steady-state error. 

A normal acceleration command (n^) equal to 1 gee is applied to the 

pitch control law which uses measurements of missile body pitch angular 
rate (q) and pitch normal acceleration (n z ) to determine the required 
actuator command ($ pc ). The actuator is modeled as a first order lag at 

188.4 rad/sec (30 HZ) and is shown in Figure 2.6. 
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Figure 2.3 Aerodynamic Sign Convention. 
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Figure 2.4 Linear Pitch Channel Dynamic Model. 
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Figure 2.6 Uncoupled Pitch Channel Autopilot. 



3. Transfer Functions o f Aerodunamlc Model 

Uncoupled pitch aerodynamic transfer functions are 

Es 

+ 1 

q k(AE-BC) k(AE-BC) 

— = , deg/sec/deg ( II.C.3-1) 

Sp -C s 2 A k s 






AE-BC 



-C 



+ 1 



-C -C 



Bs 2 



AE-BC 



+ 1 



s 2 A k s 

— + + i 



-C -C 



g's/deg 



(II.C.3-2) 



q 5 



A = 



C Na , 8 = 



W 



q S 



Cns p , C = 



(57.3) q S d 






W 



(57.3) q 5 d 



l yy 



1845 



E = 



Iyy 



p , ^ = 



Where: 



V=Ma=3.95 x 968.47= 3825.46 ft/sec 

S= 71 ft 2 

d= 2 ft 

q= 1650 lb/ft 2 

Ixx= 804 slug-ft 2 

W= 2525 lbs 

Cm a = -0.06 
8 mS p = "0.08 

Cfyj =0.15 
C N8 p = 004 
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As referred to in Appendix A and Table I. 

Then 

A =0.30793 
B =0.08211 
C = -44.33162 
E = -59.10883 
k = 0.48229 

And substituting the aforementioned values into equations (II.C.3-1) 
and (II.C.3-2) they become: 



q -1.3361 s- 0.159 

= (1/sec) (II. C. 3-3) 

Sp 22.545 x 1 0-3 s2 + 3.3634 x 1 0*3 s + 1 



n z 1 06.074 x 10-3s2- 18.82 

= (g's/rad) (II.C.3-4) 

$p 22.545 x 1 0-3 s2 + 3.3634 x 1 0"3 s + 1 



Equations (II.C.3-3) and (II.C.3-4) are the pitch aerodynamic transfer 
functions for the pitch angular rate (q) about the Y B and the achieved 
maneuver acceleration in the Z B direction. 

4. Desi g n A p proach and Analusis of the Continuous Open Loop S y stem 
A fixed flight condition at 60 KFT altitude and M=3.95 was selected 
for these preliminary performance studies of circular airframe. This flight 
condition provides a sufficiently low dynamic pressure, so that missile 
maneuvers will result in large enough angles-of-attack, to exercise side 
slip control. The above fixed condition was selected, in order to reduce the 
complexities of CBTT autopilots. 
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TABLE I 



Linearized Aerodynamic Derivatives (M=3.95 ) 
Circular Airframe 





a=0° 


8=10° 


a=20° 




-0.065 


-0.082 


-0.111 




-0.025 


-0.019 


-0.003 




0 


-0.009 


-0.020 


Y 


0.021 


0.022 


0.028 


CnS Y 


-0.050 


-0.053 


-0.062 


C 18Y 


0 


-0.016 


-0.038 


Cysr 


0 


-0.009 


-0.022 


R 


0 


0.018 


0.044 


C 18R 


0.031 


0.035 


0.044 




0.15 


0.17 


0.22 


Cnsp 


0.04 


0.04 


0.05 


^ma 


-0.060 


-0.065 


-0.1 18 


QnS P 


-0.080 


-0.095 


-0.1 15 



Reference Center of Gravity (C.G) at 0.6 1 
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The analysis of the uncoupled pitch channel autopilot is based, on the 
time responses, and Bode plots of maneuver plane acceleration, body angular 
rates and tail incidence angles, obtained by OPTSYS control program of NPS. 
For that reason the differential equations of the uncoupled pitch channel 
autopilot must be put In state-space form X=FX+GU, which Is obtained as 
follows: 

a. Equations of Control Law 

(1) Acceleration Filter Equation 

1 

X J = n 2 (II.C.4-1) 

s 

+ 1 

150 

or 

X 1 S+ 150X 1 = 150n 2 (II.C.4-2) 

and utilizing inverse Laplace Transformation it becomes: 

X, = -150 X 1 + 150 nz (II.C.4-3) 

(2) Acceleration Compensator Equation 



Y 



s + 0.143 

(-0.0387) 

0.143 



(n^-X,) s (s + 5) 



(II.C.4-4) 



5 

or by minor manipulation 

Y -0.27063 s -0.0387 

= (II.C.4-5) 

(n Z0 -X,) 0.2s 2 +s 

and utilizing inverse Laplace Transformation it becomes: 
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Y+5Y=- 1 .353 1 (ri^-Xj )-0. 1 9Z5(r\ zo -X ] ) (II.C.4-6) 

Using state-space representation of a system, in which the 
forcing function involves derivatives terms (Ref. 5: pp 675-678] one yields: 
* 2 = X 3 - 1 .353 1 (n^-X, ) (II. C. 4-7) 

X 3 = -5X 3 +6.572(n 2C -X 1 ) (II.C.4-8) 

Y = X2 (II.C.4-9) 

(3) Rate Compensator Eouation 



(- 15.6) 



s + 6.55 



6.55 



$pc - 



S + 0.143 



( X 2 - X 4 ) 



(II.C.4- 10) 



0.143 



-2.3816 s- 15.6 

6po = ( X 2 - X 4 ) (II. C. 4- 1 1) 

6.993 s + 1 



and using inverse Laplace Transformation it yields: 

Sp C =-0. 1 435^-340.58 x 10-3(X 2 -X 4 )-2.2308(X 2 -X 4 ) (II.C.4-12) 

b. Actuator Equation 
1 

$ p = V (II. C. 4- 1 3) 

s + 188.4 

188.4 

and using inverse Laplace Transformation it yields: 

5p = - 1 88.4 5p+ 1 88.4 5po (II.C.4- 1 4) 
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c. Equations of Aerodynamic Model 

From aerodynamic transfer function, Eq. (II.C.3-3) 

q -1.3361 S- 0.159 

= , or 

Sp 22.545 x 1 0 -s s 2 + 3.3634 x 1 0 ‘ 3 s + I 

22.545 x 1 0-3 q s 2 + 3.3634 x I0'3qs + q = -1.3361 SpS-0.159Sp 

(II.C.4-15) 

and using inverse Laplace Transformation it yields: 
q + 1 49 . 1 86 1 0-3 q + 44.3557 q =-59.2637 Sp -7.0526 Sp (IIX.4-16) 

Using state-space representation of a system, in which the forcing 
function involves derivatives terms [Ref. 5: pp 675-6781 one obtain: 

X 4 = X 5 -59.2637 Sp (II.C.4- 1 7) 

*5 =- 1 49. 1 86 x 1 0 ' 3 X 5 -44.3557 X4 ♦ 1 .7887 S p (II.C.4- 1 8) 

X4 = q (II.C.4-19) 

From aerodynamic transfer function Eq.(II.C.3-4) 

n z 1 06.074 x 10*3 S 2 - |8.82 

= , or 

Sp 22.545 x 1 0-3 s 2 + 3.3634 x 1 0-3 s + 1 



22.545 x 1 0-3 s 2 rij + 3.3634 x 1 0 ' 3 s rij + n z = 1 06.074 xIO " 3 s 2 5 p -1 8.82S p 

(II.C.4-20) 

or ’lij + 1 49. 1 86 x 1 0*3 tij + 44.3557 n z = 4.705 Xp - 834.775 Sp 



and proceeding as above it yields: 
Xs = X 7 - 701.92 x 1 0 _ 3 s p 



(II.C.4-21) 

(II.C.4-22) 
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X 7 = - 1 49. 1 86 x 1 0-5 X 7 - 44.3557 X 6 - 1 043.366,, (II.C.4-23) 

n* z X^ + 4.7056^ (II.C.4-24) 



Utilizing state-space representation, all the aforementioned 

equations can be modeled in a tenth order system as follows: 

X,z -150 X, + 1 50 Xj + 705.75 Sp 

X 2 z X 3 - 1.3531 rij e + 1.3531 X, 

X 3 z -5 X 3 + 6.572 rij C - 6.572 X, 

*4 = X 5 - 59.2637 6p 

*5 z - 1 49. 1 86 x 1 0'* X 5 -44.3557 X* + 1 .7887 6 p 
Xg = X 7 - 701.92 x 10-3 Sp 
X 7 z - 1 49. 1 86 x 1 0-5 X 7 - 44.3557 X^ - 1 043.36 8p 
6p = -188.4 6^ + 188.4 6^ 

6^ z -0. 1 43 6po -2.2308 X 2 + 2.2308 X 4 -340.58 x 1 0-5 X 3 + 

+ 460.839 x 10-5 n M - 460.839 x 1 0“5 X, + 340,58 x 10-5 X 5 - 
- 20. 1 84 Sp (II.C.4-33) 

lij = -( X 7 - 887. 1 24 Sp + 886.422 Spp ) (II.C.4-34) 



(II.C.4-25) 

(II.C.4-26) 

(II.C.4-27) 

(II.C.4-28) 

(II.C.4-29) 

(II.C.4-30) 

(II.C.4-31) 

(II.C.4-32) 



A continuous, linear, constant coefficient system of differential 
equations, as above, can always be expressed as a set of first-order matrix 
differential equations X = FX+GU where U is the control input to the 
system. The output can be expressed as a linear combination of the state, X, 
and the input as yz HX+JU 



The state method of representing a dynamic system is very useful 
because it standardizes the information required into three matrices, F, G, 
and H, no matter how complicated the system is. 

In our case the matrices are: 
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X, 




X. 




0 


x 2 




X2 




-1.3531 


X 3 




X 3 




6.572 


*4 




X4 




0 


*5 


x= 


*5 


G= 


0 


*6 




* 




0 


X 7 




X7 




0 










0 


«pc 




Spc 




0.4608 


n z 

- 








0 

_ 



0 = ^0 



F= 



-150 

1.353 

-6.572 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 



0 

1 

-5 

0 

0 

0 

0 

0 



0 

0 

0 

0 



0 

0 

0 

1 



-44.36 -0.1492 



0 

0 

0 



0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

1 



0 

0 

0 

-59.26 

1.789 

-0.7019 



1-0.4608 -2.231 -0.3406 2.231 
0 0 0 0 



0 

0 -44.36 -0.1492 -1043 

0 0 0 -188.4 

0.3406 0 0 -20.18 

0 0 -1 887.1 



0 

0 

0 

0 

0 

0 

0 

188.4 

-0.143 

-886.4 



150 

0 

0 

0 

0 

0 

0 

0 

0 

0 



The state-variables of interest are: 

X 4 = q : pitch angular rate 

n z : pitch normal acceleration 

Sp : pitch tail incidence 

Utilizing OPTSYS the tenth order system was analyzed in terms of 
time and frequency responses, using an input step function representing the 
“1 gee command" and trim angle-of-attack a*=0, and the corresponding plots 

were obtained. 
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Figure 2.7 shows the pitch normal acceleration (n z ) response which 
has a 0.5 seconds time constant, steady-state error 0.005 and 2.0 % 
overshoot, and therefore is according to the requirements mentioned in 
section II.A. 

Figures 2.8 and 2.9 show the required body angular rate (q) and 
control surface deflection (fy), to achieve the acceleration response. These 

responses matched those presented in [Ref. 31. 

Figures 2.10 through 2.15 show the frequency responses of normal 
acceleration (r^), angular rate (q), and tail incidence (5 p ). The relative 

stability margins are shown in Table II. Since the gain and phase margins, as 
shown in Table II are positive, it ensures that the closed loop (controlled) 



system will be stable. 



Figures 2.16 through 2.18 show the Pole-Zero Maps of n z , 


q and 5, The 


poles locations are found to be at: 




s, =-159.7 + j 18.99 


(Acceleration Filter) 


(II.C.4-35) 


S 2 =-159.7- j 18.99 


(Acceleration Compensator) 


(II.C.4-36) 


s 3 = -8.290 + j 8.059 




(II.C.4-37) 


s 4 = -8.290 - j 8.059 


(Angular Rate) 


(II.C.4-38) 


S 5 =-3.759 + j 2.515 




(II.C.4-39) 


s$ =-3.759 - j 2.515 




(II.C.4-40) 


s 7 =-0.14390000 




(II.C.4-41) 


Sq =-0.000004251 


(Actuator) 


(II.C.4-42) 


% = -0.07459 + j 6.66 


(Rate Compensator) 


(II.C.4-43) 


S 10 = -0.07459 - j 6.66 


(Normal Acceleration) 


(II.C.4-44) 



sinse all poles are in the left-half of the s-plane the system is stable. 
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PITCH NORMAL ACCELERATION, NZ 
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Figure 2.7 Pitch Normal Acceleration vs Time; Uncoupled Pitch Channel 
Continuous Open Loop System. 
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PITCH ANGULAR RATE.Q 

0 . 16 - 0.17 - 0.16 - 0.15 - 0.14 - 0.13 - 0.12 - 0.11 - 0.10 - 0.09 - 0.00 - 0.07 - 0.00 - 0.03 - 0.04 - 0.03 - 0.02 - 0.01 0.0 




Figure 2.8 Pitch Angular Rote vs Time; Uncoupled Pitch Channel; 
Continuous Open Loop System. 



43 




PITCH TAIL INCIDENCE, <5P 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 0.045 0.050 0.055 




Figure 2.9 Pitch Tall Incidence vs Time; Ucoupled Pitch Channel; 
Continuous Open Loop System. 
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MAGNITUDE - DECIBELS 
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Figure 2.10 Pitch Normal Acceleration; Uncoupled Pitch Channel; 
Continuous Open Loop System; Gain vs Frequency. 



45 




PHASE - DEGREES 

- 269.0 - 247.0 - 225.0 - 203.0 - 161.0 - 159.0 - 137.0 - 115.0 - 93.0 - 71.0 




Figure 2.1 1 Pitch Normal Acceleration; Uncoupled Pitch Channel; 

Continuous Open Loop System; Phase vs Frequency. 
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Figure 2.12 Pitch Angular Rote; Uncoupled Pitch Channel; 

Continuous Open Loop System; Gain vs Frequency. 
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Figure 2.13 Pitch Angular Rote; Uncoupled Pitch Channel; 

Continuous Open Loop System; Phase vs Frequency. 
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MAGNITUDE - DECIBELS 

— 48.0 - 43.0 - 38.0 - 33.0 - 28.0 - 23.0 - 18.0 - 13.0 - 8.0 -3 0 




Figure 2.14 Pitch Tail Incidence; Uncoupled Pitch Channel; 

Continuous Open Loop System; Gain vs Frequency. 
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PHASE - DEGREES 

20.0 - 3.0 14.0 31.0 43. 0 65.0 02. 0 99.0 116.0 133.0 




Figure 2.15 Pilch Toll Incidence; Uncoupled Pitch Channel; 

Continuous Open Loop System; Phose vs Frequency. 
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TABLE II 



Stabilitu Margins of Uncoupled Pitch Channel 

q $ p nz 

Phase Crossover Frequency (rad/sec): 56.6896 57.2 1 5 1 

Gain Margin (DB): 29.0988 7.5450 

Gain Crossover Frequency (rad/sec): 5217.1016 375.9622 

Phase Margin (Degrees): 87.9482 62.9508 
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Figure 2.16 Pitch Normal Acceleration; Uncoupled Pitch Channel; 
Continuous Open Loop System; Pole-Zero Map. 
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Figure 2.17 Pitch Angular Rote; Uncoupled Pitch Channel; 

Continuous Open Loop System; Pole-Zero Map. 
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Figure 2.18 Pitch Toil Incidence; Uncoupled Pitch Chonnel; 

Continuous Open Loop System; Pole-Zero Map. 
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The zeros of pitch normal acceleration are found to he at: 
s, = -150.00 

52 = 13.32 

53 =-13.32 

54 = - 6.550 

S3 =-0.000008756 
% = 0.000004509 

S7 = - 0.1430 
Sg = - 0.07459 ♦ j 6.66 
s, = - 0.07459 - j 6.66 

The zeros of pitch angular rate are found to be at: 
s, =-150.00 

52 =-0.119 

53 = 0.00 1 734 

54 = - 6.550 

Sg =- 0.000867 ♦ ] 0.001516 
s 6 = - 0.000867 - j 0.001516 
s 7 =-0.1430 
Sg = - 0.07459 + ] 6.66 
s, = - 0.07459 - j 6.66 

The zeres of pitch tall Incidence are found to be at: 

S( =-150.00 

52 =-0.000001322 

53 = - 6.550 

s 4 = 0.000001189 



(ll.C.4-45) 

(Il.C.4-46) 

(1I.C.4-47) 

(II.C.4-48) 

(1I.C.4-49) 

(II.C.4-50) 

(II.C.4-51) 

(1I.C.4-52) 

(ll.C.4-53) 

(II.C.4-54) 
(II.C.4-55) 
(II.C.4-56) 
(II.C.4-57) 
(II.C.4-58) 
(Il.C.4-59) 
(1I.C.4-60) 
(ll.C.4-6 1 ) 
(Il.C.4-62) 

(II.C.4-63) 

(II.C.4-64) 

(II.C.4-65) 

(II.C.4-66) 
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s 6 = - 0.07459 + j 6.66 
s 7 = - 0.07459 - j 6.66 
Sg = - 0.07459 + j 6.66 
S9 = - 0.07459 - j 6.66 



S5 =-0.1430 



(II. C. 4-67) 
(II.C.4-68) 
(II.C.4-69) 
(II.C.4-70) 
(II.C.4-71) 



5. Design Approach and Analysis of the Discrete Open Loop System 

The control of physical systems with a digital computer is becoming 
more and more common. Many new digital control applications are being 
stimulated by microprocessor technology. Among the advantages of digital 
logic for control are the increased flexibility of the control programs and 
the decision-making or logic capability of digital systems which can be 
shared with the control function to meet other system requirements. 

In order to compute the state X(t) by use of a digital computer, the 
continuous-time state equation X = FX+GU must be converted into the 
discrete-time equation X = AX+BU , assuming that the input vector U 
changes only at the equally spaced sampling instants T. 

Utilizing the Analog-to-Digital conversion of ORACLS the continuous 
time matrices F and G were converted to the discrete time matrices A and 
B using sample time T=0.0125 sec (Frequency=80 cps). 

The computer listing of the above conversion is shown in Appendix 8. 

The TRANFUNC option of ORACLS was used to analyze the dicrete open 
loop transfer function and time responses and pole-zero maps plots were 
obtained as for the continuous system. 

Figures 2.19 shows the discrete time response of pitch normal 
acceleration (r^) which has a 0.5 seconds time constant, steady-state error 
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PITCH NORMAL ACCELERATION, NZ 

-O.l 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.0 0.9 




Figure 2.19 Pitch Normal Acceleration vs Time; Uncoupled Pitch Channel; 
Discrete Open Loop System. 
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about 0.005 and 2.0 % overshoot, and therefore is according to the 
requirements mentlonted in section II.A. 

Figures 2.20 and 2.21 show the discrete time responses of pitch 
angular rate (q), and pitch tail incidence (6 p ). These responses matched 

those presented in [Ref. 3]. An input step function representing "1 gee” was 
used in all cases. 

Figures 2.22 through 2.24 show the discrete Pole-Zero Maps of n 2 , q, 
and S p The poles are found to be at: 



Zy =0.1320 + j 0.03194 


(Acceleration Filter) 


(II.C.5-1) 


22 = 0.1320 - j 0.03194 


(Acceleration Compensator) 


(II.C.5-2) 


z 3 = 0.8970 + j 0.09067 




(II.C.5-3) 


Z 4 = 0.8970 - j 0.09067 


(Angular Rate) 


(II.C.5-4) 


z 5 = 0.9536 + j 0.02999 




(II.C.5-5) 


Z6 = 0.9536 - j 0.02999 




(II.C.5-6) 


Z 7 = 0.9982 




(II.C.5-7) 


z@ = 1.0000 


(Actuator) 


(II.C.5-8) 


Zg = 0.9956 + j 0.08307 


(Rate Compensator) 


(II.C.5-9) 


Z 10 = 0.9956 - j 0.08307 


(Normal Acceleration) 


(II.C.5-10) 


Since all the poles are within the unit circle the system is stable. 
The zeros of the discrete pitch normal acceleration are found to be at 


Z\ =-0.4663 




(II.C.5-1 1) 


Z 2 = 0.1532 




(II.C.5-1 2) 


Z 3 = 0.8466 




(II.C.5-1 3) 


Z 4 = 1.1810 




(II.C.5-1 4) 


Z 5 = 0.9214 




(II.C.5-1 5) 
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PITCH ANGULAR RATE.Q 

-0.I0-0.J7 -0.16 -0.13 -0.14 -0.13 -0.12 -0.11 -0.10 -0.09 -0.00 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0.00 




Figure 2.20 Pitch Angular Rate vs Time; Uncoupled Pitch Channel; 
Discrete Open Loop System. 
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PITCH TAIL INCIDENCE, <5P 

0.000 0.005 0.010 0 015 0.020 0 025 0.030 0.035 0.040 0.045 0.050 0.055 




Figure 2.21 Pitch Tail Incidence vs Time; Ucoupled Pitch Channel; 
Discrete Open Loop System. 
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o 




Figure 2.22 Pitch Normal Acceleration; Uncoupled Pitch Channel- 
Discrete Open Loop System; Pole-Zero Map. 
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o 




Figure 2.23 Pitch Angular Rate; Uncoupled Pitch Channel; 

Discrete Open Loop System; Pole-Zero Map. 



62 



0.0 - 0.6 - 0.4 - 0.2 0.0 0.2 0.4 0.6 0.0 



o 




Figure 2.24 Pilch Tall Incidence; Uncoupled Pitch Channel; 
Discrete Open Loop System; Pole-Zero Map. 



63 



26 = 


0.9982 


(II.C.5-16) 


27 = 


1.0000 


(ll.C.5-17) 


28 = 


0.9956 ♦ j 0.08307 


(II.C.5-18) 


29 = 


0.9956 - j 0.08307 


(II.C.5-19) 



The zeros of the discrete pitch angular rate are found to be at: 



2l =" 


2.2750 


(II.C.5-20) 


22 = 


0.1534 


(II.C.5-21) 


2 3 =" 


0.1386 


(II.C.5-22) 


24 = 


0.9985 


(II.C.5-23) 


2 5 = 


0.9214 


(II.C.5-24) 


26 = 


0.9982 


(II.C.5-25) 


2? = 


1.0000 


(II.C.5-26) 


28 = 


0.9956 + j 0.08307 


(II.C.5-27) 


29 = 


0.9956 - j 0.08307 


(II.C.5-28) 


zeros of the discrete pitch tail 


incidence are found to be at: 


21 =- 


0.4658 


(II.C.5-29) 


2 2 = 


0.1532 


(ll.C.5-30) 


23 = 


0.9956 - j 0.08307 


(II.C.5-31) 


24 = 


0.9956 + j 0.08307 


(II.C.5-32) 


2 5 = 


0.9214 


(II.C.5-33) 


26 = 


0.9982 


(II.C.5-34) 


2? = 


1.0000 


(II.C.5-35) 


2e = 


0.9956 + j 0.08307 


(II.C.5-36) 


29 = 


0.9956 - j 0.08307 


(II.C.5-37) 
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6. Comparison of Continuous and Discrete Sustems 

In this chapter the classical design of the continuous uncoupled pitch 
channel autopilot, as was developed in [Ref. 3], was reviewed, and the 
corresponding discrete channel was obtained using state-space 
representation of the system and the Analog-to-Digital Conversion option 
of ORACLS. 

A comparison of the time responses in the s-domain (Figures 2.7 
through 2.9) with those in the z-domain (Figures 2.19 through 2.21) shows 
that the two systems behave the same way. In addition both systems are 
stable as it can be seen from the pole-zero maps of the continuous system 
(Figures 2.13 through 2.15) , and the discrete system (Figures 2.22 through 
2.24) , since their poles are maping according to relation 2 =& T . 

The similarity of the continuous and discrete systems allow us from 
now on to work on the discrete domain without any loss of accurasy, and 
using the computer to make things much easier. 
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III. MODERN CONTROL DESIGN AND ANALYSIS OF LONGITUDINAL 



UNCOUPLED PITCH CHANNEL AUTOPILOT 



A. GENERAL 

The design technique for the uncoupled pitch channel, as developed in 
[Ref. 31, was classical using a combination of frequency response and root- 
locus techniques. 

In classical control theory, only the input, output, and error signals are 
considered important; the analysis and design of control systems are 
carried out using transfer functions, together with a variety of graphical 
techniques as root-locus plots and Nyquist plots. The unique characteristic 
of classical control theory is that it is based on the input-output relation of 
the system, or the transfer function. 

The main disadvantage of classical control theory is that, generally 
speaking, it is applicable only to linear time-invariant systems having a 
single input and a single output. It is powerless for time-varying systems, 
non-linear systems (except simple ones), and multiple-input-multiple- 
output systems. Thus classical techniques (the root-locus and frequency 
response methods) do not apply to the design of optimal and adaptive control 
systems, which are mostly time-varying and/or linear. 

The modem trend in engineering systems is toward greater complexity, 
due mainly to the requirements of complex task and good accuracy. Complex 
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systems may have multiple inputs and multiple outputs and may be time- 
varying. 

Because of the necessity of meeting increasingly stringent requirements 
on the performance of control systems, the increase in system 
complexity, and easy access to large-scale computers, modern control 
theory, which is a new approach to the analysis and design of complex 
control systems, has been developed since around 1960. This new approach 
is based on the concept of state-space representation. 

Modern control theory is contrasted with classical control theory in that 
the former is applicable to multiple-input-multiple-output systems, which 
may be linear or non-linear, time-invariant or time-varying, while the 
latter is applicable only to linear time-invariant single-input-single-output 
systems. Also, modern control theory is essentially a time-domain 
approach, while classical control theory is a complex frequency-domain 
approach. 

System design in classical control theory is based on trial -and-error 
procedures which, in general, will not yield optimal control systems. 
System design in modern control theory, on the other hand, enables the 
engineer to design optimal control systems with respect to given 
performance indexes. In addition, design in modern control theory can be 
carried out for a class of inputs, instead of a specific input function, such 
as the impulse function, step function, or sinusoidal function. Also, modern 
control theory enables the engineer to include initial conditions in the 
design. 
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The equations of the uncoupled pitch channel autopilot have already been 
In state-space form X = AX + BU ,1n the previous chapter , and therefore 
we proceed In the modern control design of the discrete state-feedback 
autopilot. 

B. STATE FEEDBACK DESIGNED AUTOPILOT 

1. Design A p proach and Analusls of Control Law 

One of the attractive features of state-space design methods Is that 
the procedure consists of two independent steps. One step assumes that we 
have all the states at our disposal for feedback purposes. In general, of 
course, this would be unrealistic since a practical engineer would not, as a 
rule, find it necessary to purchase this large number of sensors, especially 
since he knows that he would not need them, using classical design methods. 
The assumption that all states are available merely allows us to proceed 
with the first step, namely, the control law. The remaining step is to design 
an "estimator" which estimates the entire state vector, given measurements 
of the portion of the state provided by y = HX + JU. 

The final control algorithm will consist of the control law and the 
estimator combined where the control-law calculations are based on the 
estimated states rather than the actual states. This substitution is 
reasonable and the combined control law and estimator will give closed- 
loop dynamic characteristics which are unchanged from those assumed in 
designing the control law and the estimator separately. The dynamic system 
we obtain from the combined control law and the estimator is called the 
control ler.The first step is to get a good control law. 
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The control law is simply the feedback of a linear combination of all 
the states, that is. 



The characteristic equation of the controlled (closed loop) system is 



The control law design then consists of picking the elements of K so 
that the roots of equation III.B.1-2 are In desirable locations. 

A program logic for application of Ackermann's formula to compute 
the control law is given in Appendix C. The program reads in the sample time 
T, the discrete matrices A , B and the desired pole locations in the s-plane 
and gives as output the control law matrix K. 

The program used in this work was written in WATFIV computer 
language and is given in Appendix D. 

In section III.C.5 utilizing Analog-to-Digital conversion the discrete 
uncoupled pitch channel matrices A and B were obtained. Also the analysis 
of the continuous open loop system, obtained from OPTSYS, showed the 
eigenvalues of the classical design to be at: 



U = -KX 



(III. B. 1-1) 



det[ zl-A+BK ] = 0 



(III.B.1-2) 



s t =-159.724 + j 18.9921 
s 2 =-159.724- j 18.9921 

5 3 = -8.29048 + j 8.05932 

5 4 = -8.29048 - j 8.05932 

5 5 =-3.75925 + j 2.51463 
Sg =-3.75925- j 2:51463 
S 7 =-0.14393000 



Sg = -0.00000425 



(III.B.1-3) 
(III.B.1-4) 
(III.B.1-5) 
(III.B.1-6) 
(III.B.1-7) 
(III.B.1-8) 
(III. B. 1-9) 
(III.B.1-10) 
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s 9 = -0.074593 + J 6.65959 (III.B. 1-11) 

s 10 = -0.074593 - J 6.65959 (III.B. 1-12) 

Utilizing Ackermann's formula program (Appendix D) and using 
sample time T= 0.0125 sec (80 cps), the discrete matrices A,B (Appendix B) 
and desired poles locations in s-plane the continuous open loop eigenvalues 
(equations III.B. 1-3 through III.B. 1-12) the control law was found to be: 

K = 1-0.0188 -0.0279 -0.0042 0.0290 0.0041 

0.0 0.0001 1.7732 -2.0162 0.0131] 

(III.B. 1-13) 

Utilizing Transfer Function Analysis from TRANFUNC option of 
ORACLS, and using the discrete matrices A,B and K the state-feedback 
autopilot was designed and analysed according to step responses and Pole- 
Zero Maps. 

2. Performance of Sustem 

A comparison of the discrete open loop eigenvalues (detl zl-A 1=0), 
obtained using the Compensator Transfer Function Analysis from TRANFUNC 
option of ORACLS: 



z 1 = 0.131993 + ) 0.0319377 


(III.B.2-1) 


Z 2 = 0.131993 - j 0.0319377 


(HI.B.2-2) 


Z 3 = 0.896987 + J 0.0906708 


(III.B.2-3) 


Z 4 = 0.896987 - J 0.0906708 


(III.B.2-4) 


Z 5 = 0.953625 + j 0.0299850 


(III.B.2-5) 


Z 6 = 0.953625 - ] 0.0299850 


(III.B.2-6) 


z 7 =0.998202 


(III.B.2-7) 


Zg = 1.000000 


(III.B.2-8) 
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Zg = 0.995608 + j 0.0830713 


(III.B.2-9) 


z 10 = 0.995608 - j 0.0830713 


(III.B.2-10) 


with the closed loop regulator eigenvalues 


(detl zl-A+BX ]=0) 


obtained using the Transfer Function Analysis from 
ORACLS: 


TRANFUNC option of 


Z\ =0.134462 + ] 0.0320974 


(III.B.2-1 1) 


Z 2 =0.134462 -] 0.0320974 


(III.B.2-1 2) 


Z 3 = 0.896990 + ] 0.0906457 


(III.B.2-1 3) 


Z 4 = 0.896990 - ] 0.0906457 


(III.B.2-1 4) 


Z 5 = 0.953624 + ] 0.0299994 


(III.B.2-1 5) 


z 6 = 0.953624 - j 0.0299994 


(III.B.2-1 6) 


z 7 = 0.998202 


(III.B.2-17) 


Z 8 = 1.000000 


(III.B.2-1 8) 


Z 9 = 0.995608 + j 0.0830713 


(III.B.2-1 9) 


z 10 = 0.995608 -] 0.0830713 


(III.B.2-20) 



shows no significant differences. 

Figures 3.1 through 3.3 show the discrete closed loop time responses 
of pitch normal acceleration (n z ), pitch angular rate (q), and pitch tail 
incidence (6 p ). An input step function representing "1 gee" was used. A 
comparison of the responses mentioned above, with the corresponding time 
responses of the discrete open loop system figures 2.19 through 2.21 shows 
that they are identical. 

Figures 3.4 through 3.6 show the Pole-Zero Map plots of n z , q, and Sp. 
All the poles, given by equations III.B.2-1 1 through III.B.2-20, are inside the 
unit circle and the system is stable. There are no real zeros in the system. 
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PITCH NORMAL ACCELERATION, NZ 

0.0 0.1 0.2 0.3 0.4 0.3 0.0 0.7 0.8 0.9 




Figure 3.1 Pilch Normal Acceleration vs Time; Uncoupled Pitch Channel; 
Discrete State-Feedback Autopilot. 



72 



PITCH ANGULAR RATE.Q 

- 0 . 18 - 0.17 - 0.18 - 0.15 - 0.14 - 0.13 - 0.12 - 0.11 - 0.10 - 0.09 - 0.00 - 0.07 - 0.08 - 0.05 - 0.04 - 0.03 - 0.02 - 0.01 0.00 




Figure 3.2 Pitch Angular Rate vs Time; Uncoupled Pitch Channel; 
Discrete State-Feedback Autopilot. 
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PITCH TAIL INCIDENCE, <5P 

0.000 0.005 0.010 0 015 0.020 0.025 0.030 0.035 0.040 0.045 0.050 0.055 




Figure 3.3 Pitch Tall Incidence vs Time; Uncoupled Pitch Channel; 
Discrete State-Feedback Autopilot. 
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Figure 3.4 Pitch Normal Acceleration ; Uncoupled Pitch Channel; 
Discrete State-Feedback Autopilot; Pole-Zero Map. 
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Figure 3.5 Pilch Angular Rote; Uncoupled Pitch Channel; 

Discrete State-Feedback Autopilot; Pole-Zero Map. 
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Figure 3.6 Pitch Tail Incidence; Uncoupled Pitch Channel; 

Discrete State-Feedback Autopilot; Pole-Zero Map. 



77 



A 7V 



3. Simplified State-Feedback Designed Autopilot 

In the previous section the state-feedback autopilot was designed 
using the feedback gain vector (equation III.B.M3), and therefore there 
were ten (10) returning gain loops making the system more complex. In 
order to simplify the system, zero gains were Introdused, where those gains 
were very small in the state feedback gain vector, and in doing so returning 
gain loops were eliminated and the system became simpler. 

The resulted state-feedback gain vector Is: 

K = [-0.0188 -0.0279 0.0000 0.0290 0.0000 

0.0000 0.0000 1.7732 -2.0162 0.0131] 

(III.B.3-1) 

Following the same procedure, as in the case of State-Feedback 
Designed Autopilot, the step responses and Pole-Zero Map plots were 
obtained and analysed. 

The pole locations of the simplified state-feedback autopilot are 
found to be at: 



Z 1 = 0.134450 + j 0.0321064 


(III.B.3-2) 


z 2 = 0.134458 - j 0.0321064 


(III.B.3-3) 


z 3 = 0.896805 + j 0.0911170 


(III.B.3-4) 


z 4 = 0.896805 -] 0.0911 170 


(III.B.3-5) 


z 5 = 0.953644 + j 0.0290202 


(III.B.3-6) 


z 6 = 0.953644 - j 0.0290202 


(III.B.3-7) 


z 7 = 0.998202 


(III.B.3-8) 


z 9 = 1.000000 


(III.B.3-9) 


Z 9 = 0.995608 + j 0.0830713 


(III.B.3-10) 
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2 10 = 0.995608 - ] 0.08307 1 3 (III.B.3- 1 1 ) 

A comparison of the above pole locations with those of the state- 
feedback designed autopilot (equations III.B.2-1 1 through III.B.2-20) show a 
very small shifting of the poles of the simplified state-feedback autopilot 
to the left. 

0 

Figure 3.7 shows the time response of the pitch normal acceleration 
(rij) for the simplified state-feedback autopilot. An input step function 

representing "1 gee M was used. The response shows that the steady-state 
error is greater than in the case of state-feedback autopilot (Fig. 3.1), but 
is whithin the requirements. 

Figures 3.8 and 3.9 show the time responses of pitch angular rate (q) 
and pitch tail incidence (Sp) for the simplified state-feedback autopilot. A 

comparison with the corresponding time responses for the state-feedback 
autopilot (Fig. 3.2 and Fig. 3.3) shows no significant differences. 

Figures 3.10 through 3.12 show the Pole-Zero Map plots of n z , q, 6 p . 

All the poles, given by equations III.B.3-2 through III.B.3- 1 1, are inside the 
unit circle and the system is stable. There are no real zeros in the system. 

C. ESTIMATOR DESIGNED AUTOPILOT 

1. Desi g n Approach and Anal u sis of Discrete Closed Loop Estimator 
The control law designed in the last section III.B assumed that all 
states were available for feedback. Since, typically, not all states are 
measured, an estimator must be design in order to estimate all the states, 
given measurements of a portion of them. If the state is X, then the 
estimate is X and the idea is to let U = -LX, replacing the true states by 
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PITCH NORMAL ACCELERATION, NZ 

-O.l 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 




Figure 3.7 Pitch Normal Acceleration vs Time; Uncoupled Pitch Channel; 
Discrete Simplified State-Feedback Autopilot. 
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PITCH ANGULAR RATE.Q 

-0.ia~0.17 -0.16 -0.15 -0.14 -0.13 -0.12 -0.11 -0.10 -0.09 -0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 




Figure 3.8 Pitch Angular Rate vs Time; Uncoupled Pitch Channel; 
Discrete Simplified Stete-Feedbock Autopilot. 
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PITCH TAIL INCIDENCES 
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Figure 3.9 Pitch Tall Incidence vs Time; Uncoupled Pitch Channel; 
Discrete Simplified State-Feedback Autopilot. 
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Figure 3.10 Pilch Normal Acceleration ; Uncoupled Pitch Channel; 

Discrete Simplified State-Feedback Autopilot; Pole-Zero hap. 
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Figure 3.11 Pitch Angular Rate; Uncoupled Pitch Channel; 

Discrete Simplified State-Feedback Autopilot; Pole-Zero Map. 
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Figure 3.12 Pitch Toil Incidence; Uncoupled Pitch Channel; 

Discrete Simplified State-Feedback Autopilot; Pole-Zero Map. 
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their estimates in the control law. Figure 3.13 shows an open loop estimator 
and Figure 3.14 shows a closed loop estimator. 

The characteristic equation of the closed loop estimator is: 

det[ zl-A+LH ] = 0 (III.C.1-1) 

To select L, we take the same approach as we did when designing the 
control law K. 

If we take the transpose of A-LH, we get A T -H T L T , which is the same 
form as the system matrix A-BK of the control problem. Therefore if we 
substitute A T for A, H T for B, and L T for K, we can use the control design 
results. 

Then utilizing Ackermann’s formula program (Appendix D), and using 
sample time T= 0.0125 sec, the matrices A T for A, H T for B, where 

H T = ( 0 0 0 1 0 0 0 0 0 01 

(III.C.1-2) 

assuming that we can only measure the state X 4 = q, and taking the 

desired pole locations of the estimator to be slightly to the left from the 
original open loop eigenvalues (equations III.B.1-3 through III.B.1-12), as 
shown below: 



s 1 =-159.729 + j 18.9921 


(III.C.1-3) 


s 2 =-159.729 - j 18.9921 


(III.C.1-4) 


S 3 = -8.29548 + ] 8.05932 


(III.C.1-5) 


s 4 = -8.29548 - ] 8.05932 


(III.C.t- 6 ) 


S 5 =-3.76425 + j 2.51463 


(in.c.t-7) 


S 6 =-3.76425- j 2.51463 


(III.C.1-8) 


s 7 =-0.148930 


(in.c.i-9) 
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Figure 3.13 Open Loop Estimator. 




Figure 3.14 Closed Loop Estimator. 



Sg =-0.005004 


(III.C.1-10) 


S 9 = -0.079593 + j 6.65959 


on.c.i-i i) 


s 10 = -0.079593 - j 6.65959 


(IH.C.1-12) 


we obtain the estimator gain vector : 




L T = [ -1.8784 0.0200 -0.0972 0.0219 


0.0041 


-0.0015 0.3860 0.0720 0.0007 


-0.3370 1 




(III.C. 1-13) 


Utilizing Compensator Transfer Function Analysis 


from TRANFUNC 


option of ORACLS and using the matrices A, H and L 


the closed loop 


estimator autopilot was designed and analyzed according to step responses 


and Pole-Zero Maps. 




2. Performance of System 




The closed loop estimator eigenvalues, det[ zl-A+LH ] = 0, obtained 


from Compensator Transfer Function Analysis are: 




Z\ = 0.140601 


(III.C.2-1) 


z 2 = 0.101457 


(III.C.2-2) 


z 3 = 0.897018 + j 0.0908305 


(III.C.2-3) 


z 4 = 0.897018 - j 0.0908305 


(III.C.2-4) 


z 5 = 0.953599 + ] 0.0299821 


(III.C.2-5) 


z 6 = 0.953599 - j 0.0299821 


(III.C.2-6) 


z 7 = 0.995619 + j 0.0830597 


(III.C.2-7) 


z 9 = 0.995619 - j 0.0830597 


(III.C. 2-8) 


Zg = 0.998203 


(III.C.2-9) 


z 10 = 0.999998 


(III.C.2-10) 
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A comparison of the above eigenvalues with those of the closed loop 
regulator eigenvalues det( zl-A+BK 1=0, (equations III.B.2-11 through 
III.B.2-20), shows an acceptable difference since we estimated the entire 
state vector from the measurement of just one state, namely the pitch 
angular rate (q), and we moved the desired pole locations of the estimator 
to the left of the original open loop eigenvalues, In order to get a faster 
response. 

Figures 3.15 through 3.17 show the discrete closed loop estimator 
time responses of pitch normal acceleration (n z ), pitch angular rate (q), and 
pitch tail incidence (6 p ). An Input step function representing "1 gee" was 
used. A comparison of the responses mentioned above with the 
corresponding time responses of the state-feedback autopilot Figures 3.1 
through 3.3 shows no differences. This Is due to the fact that the same 
(zero) initial conditions was used for both cases. 

Figures 3.18 through 3.20 show the Pole-Zero Map plots of n z , q and S ? 

for both the state-feedback and the estimator designed autopilots. All the 
poles, given by equations III.C.2-1 through III.C.2-10,are Inside the unit 
circle and the system is stable. 

The zeros of the estimator's pitch normal acceleration are found to be 
at: 



Zj = 0.1053 



(III.C.2-1 1) 
(III.C.2-1 2) 
(III.C.2-1 3) 
(III.C.2-1 4) 
(III.C.2-1 5) 



z 2 = 0.8970 + j 0.09080 
Z 3 = 0.8970 - j 0.09080 
Z 4 = 0.1413 



z 5 = 0.9956 + j 0.08307 
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PITCH NORMAL ACCELERATION, NZ - ESTIMATOR 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.0 0.9 




Figure 3.15 Pitch Normal Acceleration vs Time; Uncoupled Pitch Channel; 
Discrete Estimator Autopilot. 
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PITCH ANGULAR RATE.Q - ESTIMATOR 
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Figure 3.16 Pitch Angular Rate vs Time; Uncoupled Pitch Channel; 
Discrete Estimator Autopilot. 
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FITCH TAIL INCIDENCE, <5P - ESTIMATOR 

0.000 0.005 0.010 0 015 0.020 0.025 0.030 0.035 0.040 0.045 0.050 0.055 




Figure 3.17 Pitch Tail Incidence vs Time; Uncoupled Pitch Channel; 
Discrete Estimator Autopilot. 
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Figure 3.18 Pitch Normal Acceleration ; Uncoupled Pitch Channel; Discrete 
State Feedback and Estimator Autopilots; Pole-Zero Map. 
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Figure 3.19 Pitch Angular Rate; Uncoupled Pitch Channel; Discrete State 
Feedback and Estimator Autopilots; Pole-Zero Map. 
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Figure 3.20 Pitch Tail Incidence; Uncoupled Pitch Channel; Discrete 
State Feedback and Estimator Autopilots; Pole-Zero Nap. 
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Z 6 = 0.9956 - j 0.08307 


(III.C.2-16) 


Z 7 = 0.8468 


(III.C.2-17) 


Zg = 1.0000 


(III.C.2-18) 



z, = 0.9982 (III.C.2-19) 

The zeros of the estimator's pitch angular rate are found to be at: 

Z, = 0.9956 + j 0.08307 (III.C.2-20) 



z 2 = 0.9956 - j 0.08307 


(III.C.2— 2 1 ) 


Z 3 = 0.8970 + ] 0.09081 


(III.C.2-22) 


z 4 = 0.8970 - j 0.09081 


(1I1.C.2-23) 


z 5 = 1.0000 


(III.C.2-24) 


z 6 = 0.1056 


(I1I.C.2-25) 


z 7 = 0.9982 


OII.C.2-26) 


Zq = 0.9536 + j 0.02998 


(Ill.C.2-27) 



Z 9 = 0.9536 - J 0.02998 (III.C.2-28) 

The zeros of the estimator's pitch tail incidence are found to be at: 
z, = 0.9958 + J 0.08305 (III.C.2-29) 



Z 2 = 0.9958 - ] 0.08305 


(III.C.2-30) 


z 3 = 0.9536 + j 0.02997 


(1II.C.2-31) 


z 4 = 0.9536 - ] 0.02997 


(III.C.2-32) 


Z 5 = 0.1413 


(III.C.2— 33) 


Z 6 = 1.0000 


(1I1.C.2-34) 


z 7 = 0.9982 


(III.C.2-35) 


Zg = 0.8970 + j 0.09081 


(III.C.2-36) 


Z 9 = 0.8970 -] 0.09081 


(I1I.C.2-37) 



96 



3. Simplified Estimator Designed Autopilot 

For the same reason, as In the case of simplified state-feedback 
autopilot, in order to eliminate returning loops in the design of simplified 
estimator autopilot, zero gains were introduced where the gains of the 
estimator (equation III.C.1-13) were very small. Then the resulted gain 
vector is: 

U=( -1.8784 0.0200 -0.0972 0.0219 0.0000 

0.0000 0.3860 0.0720 0.0000 -0.3370 ] 

(III.C.3-1) 

Following the same procedure, as In the case of the Estimator 
Designed Autopilot, the time responses and Pole-Zero Map plots were 
obtained and analyzed. 

The pole locations of the simplified estimator autopilot are found to 

be at: 



Z 3 = 0.896885 + j 0.0938778 
Z 4 = 0.896885 - j 0.0938778 
z 5 = 0.953882 + j 0.0288601 
z 6 = 0.953882 - j 0.0288601 
Z 7 = 0.995605 + j 0.0830655 
Zg = 0.995605 - j 0.0830655 
Z 9 = 0.998202 



Z 1 = 0.140707 
Z 2 = 0.101090 



Z 10 = 0.999985 



(III.C.3-2) 
(III. C. 3-3) 
(III. C. 3-4) 
(III. C. 3-5) 
(III. C. 3-6) 
(III.C.3-7) 
(III. C. 3-8) 
(III.C.3-9) 
(III. C. 3- 10) 
(III.C.3-1 1) 
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A comparison of the aoove pole locations with those of the estimator 
designed autopilot (equations III.C.2-1 through III.C.2-10) shows no 
significant differences. 

Figures 3.21 through 3.23 show the simplified closed loop estimator 
time responses of pitch normal acceleration (n z ), pitch angular rate (q), and 
pitch tall Incidence ($p). An Input step function representing “1 gee" was 
used. A comparison of the above responses with the corresponding time 
responses of the estimator designed autopilot (Figures 3.15 through 3.17) 
shows no significant differences. 

Figures 3.24 through 3.26 show the Pole-Zero Map plots of n z , q and S ? 

for the simplified estimator designed autopilot. All the poles, given by 
equations III.C.3-2 through III.C.3-1 1,are inside the unit circle and the 
system Is stable. There are no significant differences In the zero’s 
locations from those of the estimator autopilot (equations III.C.2-1 1 
through III.C.2-37). 
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PITCH NORMAL ACCELERATION, NZ 
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Figure 3.21 Pitch Normal Acceleration vs Time; Uncoupled Pitch Channel; 
Discrete Simplified Estimator Autopilots. 
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PITCH ANGULAR RATE.Q 
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Figure 3.22 Pitch Angular Rote vs Time; Uncoupled Pitch Channel; 
Discrete Simplified Estimator Autopilot. 
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PITCH TAIL INCIDENCE, <5P 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 0.045 0.050 0.0 




Figure 3.23 Pitch Tail Incidence vs Time; Uncoupled Pitch Channel; 
Discrete Simplified Estimator Autopilot. 
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Figure 3.24 Pitch Normal Acceleration ; Uncoupled Pitch Channel; 

Discrete Simplified Estimator Autopilot; Pole-Zero Map. 
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Figure 3.25 Pitch Angular Rate; Uncoupled Pitch Channel; 

Discrete Simplified Estimator Autopilot; Pole-Zero Map. 
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Figure 3.26 Pitch Toil Incidence; Uncoupled Pitch Channel; 

Discrete Simplified Estimator Autopilot; Pole-Zero Nap. 
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IV. DISCRETE COUPLED PITCH AND ROLL CHANNEL AUTOPILOTS 



A. GENERAL 

The control design procedures described in chapter II and III were applied 
to systems with one input and one output. The transfer-function approach in 
[Ref. 3] is fundamentally limited to single input/output systems while the 
state-space methods of chapter III were limited to single input/output to 
simplify the procedures. In fact, if one tries to apply the "pole-placement” 
approach of chapter III to a multivariable (more than one input or output) 
system, all the procedures work but the gains, K or L, are not uniquely 
determined by the resulting equations. Until a design approach is available 
which intelligently uses this extra freedom, the pole-placement concept for 
estimator designs for systems with more than one output or for controller 
designs for systems with more than one input has limited value. 

In this chapter the pole placement and robustness design program 
(POPLAR), [Ref. 8] was used as a tool in the design of multivariable systems. 
This program is designed to employ singular value analysis and the use of an 
optimization to aid in pole placement control system design for linear 
multivariable systems. Robustness of the system is also considered by 
establishing singular value levels which correspond to multiloop gain and 
phase margins determined from the universal gain phase diagram developed 
by Newsom and Mukhapadhyay at NASA Langley. 
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B. DESIGN APPROACH AND ANALYSIS OF THE DISCRETE STATE 
FEEDBACK AUTOPILOT 

The first step in any multivariable design should be an attempt to 
find an approximate model consisting of two or more single input/output 
model or else decouple the control law matrix K and the estimator law 
matrix L. This step will give better physical insight into the important 
feedback variables and may lead to a plant description which is 
substantially simpler for design purposes and yet yields no significant 
degradation from an analysis based on the full multivariable system. 

Combining the pitch and roll channel autopilot, the roll channel taken 
from [Ref. 7], a seventeenth (17* h ) order system is obtained and is shown in 
Appendix E. 

Utilizing Transfer Function Analysis from TRANFUNC option of ORACLS 
and using the combined pitch and roll matrices A, B and K the discrete 
closed loop system was obtained and analyzed according to time responses 
and Pole-Zero Map plots. 

C. PERFORMANCE OF SYSTEM 

The discrete state-feedback eigenvalues, det[ z!-A+BK]=0, of the 
coupled system, obtained from Transfer Function Analysis option of 
ORACLS, are: 



z, = 0.1 34462 ♦ j 0.0320974 


(IV.C-I) 


z 2 = 0.134462 - j 0.0320974 


(IV.C-2) 


z 3 = 0.896990 ♦ j 0.0906457 


(IV.C-3) 


z 4 = 0.896990 - j 0.0906457 


(IV.C-4) 


Z 5 = 0.953624 ♦ j 0.0299994 


(IV.C-5) 
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Z 6 = 0.953624 - j 0.0299994 


(IV.C-6) 


Z 7 = 0.998202 


(IV.C-7) 


Z 8 = 1.000000 


(IV.C-8) 


Z 9 = 0.995608 + j 0.0830713 


(IV.C-9) 


Z 10 = 0.995608 - j 0.0830713 


(IV.C-IO) 


Z\ \ - 0.0993086 


(IV.C-1 1) 


Z 12 = 0.836553 + j 0.306107 


(IV.C-1 2) 


Z 13 = 0.836553 - j 0.306107 


(IV.C-1 3) 


z 14 = 0.967788 + j 0.0334452 


(IV.C-1 4) 


Z 15 = 0.967788 - j 0.0334452 


(IV.C-1 5) 


Z 16 = 0.895730 


(IV.C-1 6) 


Z 17 = 0.938133 


(IV.C-17) 



Since all the poles are inside the unit circle the coupled system is 
stable. 

Figures 4.1 through 4.3 show the closed loop time responses of pitch 
normal acceleration (n 2 ), pitch angular rate (q) and pitch tail incidence (6 p ). 

An input step function representing "1 gee” was used. All the responses are 
identical to those obtained for the uncoupled pitch channel and therethore 
they meet the desired requirements. 

Figures 4.4 through 4.6 show the Pole-Zero hap of n 2 , q and S p 
respectively. The poles of the system are given by equations IV.C-1 through 
IV.C-17 and they are the same as the uncoupled closed loop poles (equations 
III.B.2-1 1 through III.B.2-20). There are no real zeros in the state feedback 
coupled system. 
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PITCH NORMAL ACCELERATION, NZ 
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Figure 4.1 Pitch Normal Acceleration vs Time; Coupled Pitch and Roll 
Channels; Discrete State-Feedback Autopilot. 
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PITCH ANGULAR RATE.Q 
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Figure 4.2 Pitch Angular Rote vs Time; Coupled Pitch and Roll Channels; 
Discrete State-Feedback Autopilot. 
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FITCII TAIL INCIDENCE, <5P 
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Figure 4.3 Pitch Toil Incidence vs Time; Coupled Pitch and Roll Channels; 
Discrete State-Feedback Autopilot. 
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Figure 4.4 Pitch Normal Acceleration ; Coupled Pitch and Roll Channels; 
Discrete State-Feedback Autopilot; Pole-Zero Map. 
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Figure 4.5 Pitch Angular Rate; Coupled Pitch and Roll Channels; 

Discrete State-Feedback Autopilot; Pole-Zero Map. 
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Figure 4.6 Pitch Tail Incidence; Coupled Pitch and Roll Channels; 
Discrete State-Feedback Autopilot; Pole-Zero Map. 
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A 



D. ROBUSTNESS OF SYSTEM 

With the rising interest In multivariable control theory brought on by 
increasingly complex systems the need has arisen to develop design methods 
that will allow the designer to specify system performance while at the 
same time ensuring relatively high stability margins or robustness. In the 
single-input-single-output (SISO) case the designer has had the tools to do 
these tradeoffs In the form of Nyqulst, Bode and root-locus plots. In the 
multi-lnput-multl-output (MIMO) case the classical methods are not totally 
appropriate. 

With the Increased Interest In MIMO systems numerous methods of design 
have been employed to obtain suitable system performance and robustness 
with varying degrees of success. One primary method of design is to keep 
the plant as decoupled as possible throughout the design so that each 
Individual element may be controlled Independently and designed essentially 
as a single loop system. In another procedure the multiloop system is 
modified Into a system that has diagonal elements that are much larger than 
any off-diagonal elements. This diagonally dominant system Is then in a 
form where conventional Nyquist type techniques can be employed in the 
analysis. A third common MIMO design method Is that of the Linear Quadratic 
(LQ) method. This method uses a quadratic cost functional and 
optimamization principles to allow the designer to design for various 
performance levels by adjusting the matrix weighting terms used In the 
cost function. The major difficulty with all of the above methods is that 
they are not necessarily robust. This Is especially true for cross-coupling 
terms between loops. 
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A simple Interpretation of robustness Is the ability of a system to 
tolerate design pertubatlons. These pertubatlons could be in the form of 
actuator failures, plant parameter uncertainty, unmodeled dynamic or 
nonlinear terms, or any one of many other pertubations to the nominal 
design of the system. 

The analysis of the robustness of the system was obtained using the Pole 
Placement and Robustness Design Program (POPLAR) developed In [Ref. 8]. 

Utilizing the above program and using the data file shown In Appendix F 
for the state feedback autopilot, and the data file shown in Appendix G for 
the estimator autopilot, the system found to be robust. 

Figures 4.7 and 4.8 show the Minimum Additive Output Singular Value 
(SVADMO) vs frequency for the state feedback and estimator designed 
autopilots respectively. It Is noted that for very low frequencies the values 
of SVADMO are above 0.80680 and 0.83976 for each of the designs. 

Figures 4.9 and 4.10 show the Minimum Additive Input Singular Value 
(MIN.ADD.IN.SV) vs frequency for the state feedback and estimator designed 
autopilots respectively. 

In the case of state feedback autopilot the ordered complex eigenvalues 
(Input), and the ordered computed eigenvalues (output) are: 

QPte rad . ln oyl Ordered out put 



-0.01000 + j 1.00000 
0.09928 

0.13199- j 0.03194 
0.13199 + j 0.03194 
0.83522 - ] 0.30974 



0.09931 

0.13446 - j 0.03210 
0.13446 + j 0.03210 
0.83653 - j 0.30610 
0.83653 + j 0.30610 
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SVADMO 

0.800 0.025 0.050 0.875 0.900 0.925 0.950 0.975 1.000 1.025 1.050 




Figure 4.7 SVADMO vs Frequency; Coupled Pitch and Roll Channels; 
Discrete State-Feedback Autopilot. 
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SVADMO 

0.039760.05976 0.07976 0.09976 0.91976 0.93976 0.95976 0.97976 0.99976 1.01976 




Figure 4.0 SVADMO vs Frequency; Coupled Pitch and Roll Channels 
Discrete Estimator Autopilot. 
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MIN ADD IN SV 

0.999a. 0000 1.0001 1.0002 1.0003 1.0001 1.0005 
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o 




Figure 4.9 MIN. ADD. IN. SV vs Frequency; Coupled Pitch and Roll Channels; 
Discrete State-Feedback Autopilot. 
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MIN ADD IN SV 

0.91 0.95 0.9G 0.97 0.90 0.99 1.00 1.01 




Figure 4.10 MIN.ADD.IN.SV vs Frequency; Coupled Pitch and Roll Channels; 
Discrete Estimator Autopilot. 

1 19 



0.83522 + j 0.30974 
0.89380 

0.89699 - j 0.09067 
0.89699 + j 0.09067 
0.95363 - j 0.02998 
0.95363 + j 0.02998 
0.96909 - j 0.03286 
0.96909 + j 0.03286 
0.99561 - j 0.08307 
0.99561 + ] 0.08307 
0.99820 
1.00000 



0.89571 

0.89698 - j 0.09064 
0.89698 + ] 0.09064 
0.93813 

0.95360 - ] 0.03001 
0.95360 + j 0.03001 
0.96776 - ] 0.03345 
0.96776 + j 0.03345 
0.99559 - ] 0.08307 
0.99559 + ] 0.08307 
0.99820 
1.00000 



In the case of the estimator autopilot the ordered complex eigenvalues 
(Input) are the same, as In the state feedback autopilot, and the ordered 
computed eigenvalues (output) are: 

Ordered output eigenvalues 

0.07717 

0.09855 

0.25923 

0.69508 

0.83707 + j 0.30140 
0.83707 -) 0.30140 
0.89849 
0.90064 
0.93996 
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0.96207 

0.96562 + J 0.03401 
0.96562 - j 0.03401 
0.99559 + j 0.08307 
0.99559 - j 0.08307 
0.99820 
1.00000 
1.06596 



121 



V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The goal of this thesis was to design and analyze the discrete 
longitudinal autopilots for BTT and STT missiles. 

The following are the principal conclusions based on this work. 

1. The continuous and discrete longitudinal autopilots were proved to 
have identical performances. 

2. The State-Feedback and Estimator autopilots were introduced as 
additional dynamic designs in order to implement control into the original 
system. Both performances were found to meet the desired requirements. 

3. The Simplified State-Feedback and Estimator autopilots reduced the 
complexity of the system, by eliminating returning gain loops, and without 
any significant effects on the system's performance. 

4. The performance of the Coupled Pitch and Roll channel autopilots was 
found to be satisfactory and also proved to be robust. 

B. RECOMMENDATIONS 

In order to make the whole design even simpler more returning gains 
loops, of the Simplified State-Feedback and Estimator autopilots, must be 
eliminated. A further investigation then must be conducted in order to 
examine whether the performance of the resulting design remain unchange. 
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APPENDIX A 



GEOMETRIC AND MASS PROPERTIES OF MISSILE 


CONFIGURATION 


Weight (lbs) 


2525 


Maximum Diameter (in) 


24 


Length (in) 


160 


I xx (Slug Ft*) 


40 


Iyy (Slug Ft*) 


804 


l zz (Slug Ft 2 ) 


810 


Center of Gravity Distance from Nose (in) 


100.8 


Reference Length (Ft) 


2 


Reference Area (Ft 2 ). 


71 
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OPIMAT DATA F ILc INPjI TO UKLUMV PkUGKA^: OKACLr.S P^JJRAMS 



APPENDIX B 



ANALOG TO DIGITAL CONVERSION OF UNCOUPLED PITCH CHANNEL 

AUTOPILOT 
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APPENDIX C 



PROGRAM LOGIC FOR COMPUTING CONTROL LAW K VIA ACKERMANN’S FORMULA 

1. Read In A, B, T, and Nj, the number of states. 

2. Comment: first we will read in the desired pole location in the 
plane, convert them to z-plane polynomial coefficients, and construct 
<x(A). 

3. I + Identity matrix, N f x N s 

4 ALPHA + I 

5. k + 1 

6. If k > Nj, go to step 18. 

7. Read in pole location k as a+jb. 

8. If b=0, go to step 14 

9. A t + -2 exp(aT)cos(bT) 

10. A 2 + exp(2aT) 

11. ALPHA + ALPHA X (A X A + A t A +A 2 D 

12. k + k + 2 

13. Go to step 6. 

14 At + exp(aT) 

15. ALPHA + ALPHA X (A - At X I) 

16. k <r k + 1 

17. Go to step 6. 

18. Comment: now we construct the controllability matrix. 

19. C I 
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20. E * B 

21. k + 1 

22. If k > N s , go to step 28. 

23. Comment: replace column k of C by E. 

24. C[ ; k] E 

25. k + k + 1 

26. E + A X E 

27. Go to step 22. 

28. Comment: now solve for the control law, first form e T n as the last 
row of I 

29. E * UN,; 1 

30. Solve DC = E for D. 

31. K = D X ALPHA 

32. END 
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APPENDIX D 



COMPUTER PROGRAM FOR APPLICATION OF ACKERMANN'S FORMULA 



FILE: ACKERMAN * AT r IV A1 



l J Co 

INI EGE K N S» R » I A * I Ea * I » J * N I 

REAL Ph!2C,2o ) ,GA l 20,2 J) »AI ( 20, 20), ALPHA (20,23) . *i\< 40) 

REAL A2 ,c < 20 ,40 ) ,0 ( 20, 20) , 0 120) , oMc 0 ,2 )) ,JRi 20. 2 0) 

REAL A ,0 , T , R, ATE HP (20, 20) 

REAL «1 ,TG(20,20),TI(2J,20),TS(20,20),v.<20,20),TM(20,20) 

INPUT OF N S= NUMBER uF STaTES 
T= Sample TIf’c 

PRINT , • ENTER NUMBER OF STATES.' 

REAU, NS 

P R I NT , • ENTER SAMPLE TIME • 

READ, T 

INPUT JF P nl = A MATRIX 
OU 100 1=1, NS 

DO 101 J= 1 , S 
WRI T E ( 6 ,103 ) I,J 

103 FORMAT ( 4X, ' c.N IE* P H< ' , 1 2 , ’ , • , 12 , ' ) ' ) 

READ, PHI I,J ) 

101 CONTINUE 
luo Continue 

INPUT OF G A=3 MATRIX 
DO 200 1=1 ,NS 

J = l 

WRI TE ! 6 , 105 ) 1,0 

105 FORMAT < 4X, • cNTEr, GA < • , 12 , • , ' , 12, • ) ' ) 

READ ,GA ( I , J ) 

200 CONTINUE 

BUILDING THE IDENTITY MATRIX 
AlPHA=C= 1 

DO l 1 = 1, NS 
Do 5 J= L.NS 
IF ! I .E O.J ) GC TO 300 
Alt I , J ) =0 . 

ALPHA < I , j ) = 0 . 

C < I , J) = 0. 

GL To 5 

300 A I ( I , J ) =1 . 

ALPHA! I , J ) = L . 

C ( I ,J)= l. 

5 CONTINUE 

1 CONTINUE 
R=1 

16 IF (K.GT.NS) GO TO Id 

INPUT OF DEolREO PULE., LOCATION 
WRITE ( 6,43 ) 

43 FORMAT! • ENTER REAL PART OF DESIRED POLE LOCATION.') 

READ, A 

WK I TE ( 6 ,44 ) 

44 FORMAT ( ' cNTER IMAGINARY PART OF DESIRED POLE LOCATION.') 

READ.o 

C BUILDING ALPHA (PH I) 

IF ( B . EJ.O ) GO TO 14 
A 1= -2 . * tXP ( m* T) *COS( «J*T ) 

A2= EXP ( 2.*a*T ) 

C EVALUATE PPI*PHI 

CALL VMULFF < PH ,PH , HS , NS , NS , 20 ,20 , TM ,2 0, I ER ) 

OC 2 1 = 1,. NS 
OJ o J= 1 , NS 

TO ( I , J) =Ph ( I, j ) * A 1 
T I ( I , J) =A I ( I, 0 ) = A2 
T S ( I , J ) = T M ( I , J ) •*• T U ! I * J ) ♦ T I (I,J) 

6 CONTINUE 

2 CCNTINUE 

CALL VMULFF l ALPhA , TS . NS , NS , NS ,20 , 20 ,A TEMP, 20 , IE R ) 

DO 600 1=1, NS 
00 601 J=l ,<NS 

601 ALPHA! I , J ) = A TEMP ! I , J ) 

6C0 CCNTINUE 



AC RODS' ID 
AC ROlO 2D 
<-» C P\ OC J 3 V, 
AC nO DO *» D 
AC N uOQ 50 
ACRlwOoU 
AC m)oO TO 
AC mjuO SO 
ACRCC09D 
AC ROwl 00 
ACRCC l ID 
ttCr.0wl2D 
ACnOv.1 30 
mCROOIhO 
AC /xOol 50 
->C xO , l 63 
AC .\00 1 7 J 
AC Xu 31 BO 
AC nOoI 9C 
AC xCu2 00 
AwRCu2 10 
ACR0v,220 
AC.v0C230 
AC i\0o2 40 
ACR0G250 
*xC RC 0 260 
A „ 0 ^2 7C 

^ w R 0 w. 2 GO 

ACRJJ2 30 
|\C J3C0 

ac <uo3 10 

**C<'JJ320 
AO <\ J J3 30 
AC.<00 j AO 
A C<\0 03 50 
A Z \ 0 j 3 6 0 
} C .< 0 ^ 3 7 vJ 
o 0 

a: MJO jSO 
AC KOuAOO 
AC 10 
AC < C 0 4 2 C 
AC K. J j A 30 
AC \0j4h 0 
AC 0 

AC < J • j 4 6 0 
AC \ J u4 7 0 
AC < 0 j4 8 J 
AC <0 34^0 
AC, <00500 

A- i\ uj 5 10 

AC < G05 20 
AC <0 J 5 30 
AC <00 5 AO 
AC <0 05 5 0 
ACK00560 
AC <i)0 5 70 
AC <00580 
A C <\ Q 0 i? 9 0 
AC K 0u6 00 
A C < 0 0 6 10 
A3 K0u620 
ACn0w630 
AC < 0 • j 6 4 0 
AC \JUo 50 
AC K 0 06 6 0 
AC <0l>6 70 
AC < 0 ofc oO 
AC KO 06 VO 
AC <00700 
a: <007 10 
AC <00 7 20 
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FILE: ACKERMAN <•< A T r IV 



A l 







K = n+2 


a: ?7 io 






GO TC 1c 


A3 \uu7 ^ 0 




1<* 


A 1= c XP l A a T ) 


*C<00750 






DO 3 I = l * N S 


ACnOu7o3 






DO 7 J= 1 ,NS 


AC.\Gj770 






T G t I , J) = Al ( I , J) * A l 


a: ,\C J7 30 




7 


TIII.JI =PH< I , JJ-TCl I , J ) 


AC <00 7 90 




3 


CONTINUE 


mC < C Go 00 






CALL VMULFF t ALPHA, T I , N5..NS ,N a , 20 , 00 , AT EMP.2 3, 15 R J 


A l K0c3 10 






CO 700 1=1, NS 


40 rx 00820 






DC 701 J=1,NS 


AO <00630 




701 


ALPHA t I ,JI = ATlMPII ,J) 


AC<uu3A0 




700 


CONI INUE 


AOn 00850 






K=Kf 1 


A 0 K 0 0 9 6 0 






GO TO 16 


\0 <0 06 70 




19 


DC 9 1 = 1, NS 


AuKOuSSO 






J = 1 


mC KC03 90 




9 


E(l»J) = GAll,J) 


AO <\ C C 9 G 0 






K = 1 


A 0 .\ C 0 9 10 




22 


IF IK.GT.M3 ) GO TO 28 


AC i\009d J 






DO 10 1=1, NS 


40<C j9 30 






C (I ,K) = EU , 1 ) 


AG rx0o9 AQ 




10 


CONTINUE 


A 0 f\C09*30 






N=K «■ 1 


AC /< Ow’9 oO 






CALL VMULFF (PH.E.NS, NS. 1 , 2 J , 2 J , c , 20, IEK ) 


A0K3O97G 






GO TC 22 


aO < 0C9 60 




28 


DO 11 J = 1 » N S 


AC<0u*990 




11 


El 1 , J 1 =AI | NS ,J ) 


4C<O10C0 






R =0 . 


< 0 L 0 10 






Call iu irir ic,2J,i*S,',s,ivfO,23,3f^ii, lEx) 


AC<01020 


c 






AO. <01030 






CALL VMULFF lE.J.l , NS , NS , 20 , 20 , 3f , 20 , IER ) 


AC*\JiO*tO 






CALL VMULFF IdF.^LPHA.I , NS , NS ,20 , 2 0 , J K, 20 , IE PJ 


AO <0 1 J 50 


c 






AO < vj 1 0 60 






PRINT, ' CONTROL GAIN VECTOr * 


AC <0 10 70 






WRITE 16, A5) l UKl 1, J J ,J-1,NS ) 


AO. <010 sO 




4 5 


FORMAT (BP 15.4) 


A0401U9Q 






RETURN 


AOrxOl 100 






END 


AOnOII 10 


SENTRY 




AwK 01120 
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APPENDIX E 



ANALOG TO DIGITAL CONVERSION OF COUPLED PITCH AND ROLL CHANNEL 

AUTOPILOTS 
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APPENDIX F 



ROBUSTNESS DATA FILE FOR THE DISCRETE STATE FEEDBACK DESIGNED 

AUTOPILOT 



FILE: FEEDBACK DATA AL 

1000 0 



200.0 


0 


.0 1.0 








1 




0 0 








0.4 


0 


.4 0.0 0 








0 0 0 




533 10CO 


0 






l 717 170217 170C 171 7 










1.5 795758E-C1 




3. 63 69 3Z6E -0 2 


5. 58485 92E -03 


-3. 54152866-2 


-5. 6 3 5063 4 6-3 


1.437 6905E— 03 




-o. 35:3 01 9 dE- 3 


3 . 12 1 50 380 E + Q 


-2. 74 346 96 E +0 


8 .49 ou 3 4 16-01 


0 

0 




0 

0 


0 


0 


0 


7.3 7021 2 IE- 03 




1 .0001 84 76-0 0 


1 .2145713E-02 


-1. 81108296— 4 


-2 .37079035-5 


6.4803277E-C6 




— 4 • 314B36 7 E - 5 


2 • 35 544 1 5c -02 


-2. 1657925c— 2 


9. 078 .7436-03 


0 

0 




0 

0 


0 


0 


0 


-3. 5755072E-2 




- 8. 968 35 5 6 E-t 


9.3927554E-01 


8. 7928240E-04 


1 .39377 5 5 c - 04 


-3. 1463045E-5 




2 .094o99oc-04 


- 1. 14342466-1 


1 .0513^9^-01 


-4 .4066 5 5 1 1 -2 


0 

0 




0 

0 


0 


0 


0 


6.3641914E- 04 




4.8 0 723*t 76 — 0 3 


7.3780452E-04 


Q. U 345 92 t - 0 1 


1 . 1 7240 94 c- 02 


1.03891 730E-7 




- l.J2 7o90oc-o 


-2.5772733E-1 


— *♦ . 3 7o06 7o E - 1 


3.o239644c — 04 


0 

0 




0 

0 


0 


0 


0 


-1. 2632605E-* t 




— 8. 3 10 21 3 3 E “ ^ 


- 1.35084525-4 


-5. 52425915-1 


9.94jl431c-01 


-i. 1935039E-8 




1 .3 467 7 57E-0 7 


1 . 10524796-01 


1.03735536-01 


—5 .65300o2E— 5 


0 

0 




0 

0 


0 


0 


0 


5. 0 C58503E- 05 




3 . 4 906 7 02 E “0 4 


5. 352I335E-05 


-3.423525 7 5 -4 


-5.42321946-5 


9. 9653837E- 01 




1 .2 ^ 73 6736—02 


— 4 . 3 8 3 7 9 60 6-2 


- 4 . 3103 1^96-2 


2 .239022306-5 


0 

0 




0 

0 


0 


0 


0 


1. 1 198072E-02 




8 .458971 5E-02 


1 .29826646-02 


-8. 25783366-2 


-1 . 3194h33E-2 


-5. 53288 15E-1 




9 . 94659846-01 


-4. 531317^6-0 


-7. 69 330 73 E -0 


6. 37744246-03 


0 

0 




0 

0 


0 


0 


0 


-1. 8197865E -3 




-1. 64759365“ 2 


-2.5342871 5-3 


1. 54 144375-02 


2 . 59274896-03 


-7.6937266E-7 




o • L 1 29 2 60 E — 0 6 


1 . 9o 5 32 66t -02 


8.2551326 — 0 i 


- 1 • o 1 1 7 69 i E-3 


0 

0 




0 

0 


0 


0 


0 


-2. 4201 1 18E-3 




-2.6250703E-2 


- 4 . 0 4u8 ill E-3 


2 .50737846-02 


4.1o749 756-03 


-2. 2085ol2E-6 




1 .**o5903^E-05 


-1. J1 J45586-1 


8 .52133305-01 


-3 . Gou6 4 3 76 — 3 


0 

0 




0 

0 


0 


0 


0 


8.5120361E-03 




7.7170206E-02 


1 . 1870299E -02 


-7. 45350O6E-2 


- 1 .2144650E-2 


3.464751 BE - 03 




- 1. 2502&40E-2 


4 • 6 5o34 625 -Ou 


-3. 34 J95 97E-0 


1 . Ou 7 5 o 1 OE - 00 


0 

0 




0 

0 


0 


0 


0 


0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


9.97738 78E- 01 




-2. 7 Cl 36 9 OE- 3 


1 .02458966-03 


2 .20857366-01 


-2.2392473E-1 


3.22851 5 5E — 00 




3.O7504J3E-J2 








0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


1. 249 25 12E-02 




9.99991216-01 


3. 36863396-Go 


7. 92138815-04 


-7.8251 019E-4 


1.612 06 056— 02 




3 .2 1 4o 565E-J 4 








0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


-1.32996046-3 




-2. J935729E- 1 


9.0433727E— 01 


-3. 68819056-5 


3.04511 286-05 


—9. 65 97 3 1 8E —4 




-2. o 90837 0 E-5 








0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


-9. 07977236-3 




- 1. 102403 7E-2 


3 .99972925-03 


9.2 7 6303 2E- 0 1 


1.1 630732E-02 


-1. 70C79726-1 




-2. 0028512E-3 








0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


-1. 7340942E — 4 




— 2.07178 99E-4 


7.85 645 38E -05 


1.7 37 4049E - 02 


9. 10594 31E-01 


2. 450551 OE- 01 




2.725751 JE-03 








0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


-1. 9081030E-3 




-2. 30o 1830E-3 


8.4o09805E-04 


1 .5230715t-01 


-1.5037418E— 1 


7 .469 27 4 7E - 01 




- 1. 14294856-3 








0 




0 


0 


0 


0 


0 




0 


0 


0 


0 
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FILE: FEEDBACK DATA 



A I 



-6.2875718E-2 
4. 3 32 3648E + C 1 
-7. 56O9029E-3 
-1.6449186E-2 
7.9d2 1623t-C2 
-9. 98 01 608 1 -4 
1.82856 12E-04 
•7.2440 726c -5 
-1. 7 5 75 5 15 1 -2 
3.<t31 5556E-C3 
5.4307555E- 03 
-1 .o073019E-2 
0 
0 
0 
0 
0 
0 
0 

1. 00000 00E-C0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

1. OCOOOOOE- CO 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

1.000 00 OOE- 00 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



-7,3 oC 2 h9 7E- 2 
5 • 1 J 4J 68 5 E- 02 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

2 .7 C 13 64 OE-O 3 
6 • 7 92o 9 7 1 e -0 t> 
2 .0 935 72 9E-C 1 
1 .10240675-02 
2 .07178995-04 
2 .3081 93 OE-O 3 
7 . 5o02 497E-0 2 
0 
0 
0 
0 

1 .OCOCOOOE-OO 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

1 . COCO 0 GOE-3 0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

1 . 0 COO OQOE-3 0 
0 
0 
0 
0 
0 
0 



2 .81 7ol912-Q2 



0 

0 

0 

0 

0 

0 

1 .00000005-00 
0 
0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 • CJOOO OOE -00 
0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 .00000005-00 
0 



D.ol 22435 l-0G 



0 

0 

0 

0 

0 

0 

0 

0 

0 

1 . JQ00Q00E-00 
0 
0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 . OOOOOOOc-OO 
0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



-5.5314j55E-G 



1.00003 OOfc- 



1 . 00U03 OOE - 



0 
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ooo ooo ooo ooo ooo ooo ooo ooo ooo ooo ooo ooo ooo 



FILE: F E El) d AC K OAT.; 



A L 



1 . OCQOOGCE 



- 0 . 013800 GE 
O.OUOGOOOE 
0 . GGGCOOOE 
O.OOGGOOOE 
O.OOOOOOOE 
C. OOOOOOOt 
C.OCOOOOGt 
- 1.4 990000 E 
0 G 0 0 0 G 
0 C 0 0 0 G 
-10000 
0.13199300 
0.13199300 
0.89698700 
0 .89698700 
0.95362500 
0.95362500 
0.99 82 02 00 
1.00000000 
0.99560800 
0.99560800 
0 .09 -*2 805 0 
0 .33521600 
0.83521600 
0 .96908700 
0 . ^o 9 C 8700 
C . 39379700 
0.93707500 



0 O 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

-GO 0 

0 0 

0 0 

0 0 

0 l . OCCOOOOE-J 0 

-00 -0. J279GC0E-CG 

>00 0. OCCIQOOE>00 

♦ GO U.OCCOG J0E>00 

>C0 0. OOOOOOOE+O 

>00 O.OCCOOGQE>OQ 

-GO J.OCCOOOOc>OU 

> CO 0 • OOG-t-OOOr -00 

-00 O.GC23000E-00 

OOOQOCuUOO 
COO COCOOOO 
10000 

0.03 1937 70 
-G.0JI9377G 
0.09067030 
-0.09^67030 
G . 02 9935 00 

- C. 02998500 
0 . GOOOOOOG 
C.OOJCCOUO 
0.03 3071 3C 

-C. 083 0 71 30 
v«. OOOGOO uu 
C. J>0 ) 7 56 30 

- C . 5 0 'f 7 36 0 0 
G. 03285380 

-C. 032 35330 
C. 000000 JO 
C.OOOCLOOQ 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0 . 00420006-00 
1 . 77320 OOE >00 
O.OOOGDOoc *00 

O.COOOOOOE>OG 
O.OOwOOUOE>Q 
0 . C 3 Q 0000 E +00 



0 

L . JJOOOOQc- oO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 02900 OOoOO 
-2.0162000000 
0. 00O00G0E>00 

0 . OOOOOOjE >0 
0. 0 0 GOGGLE >0 
0. 01 78000 e -0 



0 

0 

0 

0 

1. 0000 3 OOE - 00 

0 

0 

0 

0 

0 

0 

0. 0091 0 CGE >00 

o.eu i oooe>oo 

0. QG00Q00E> 00 

0 . GOO J G OOE > 00 
C . GOuO OGoE >0 
-0.0179000E-0 
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APPENDIX G 



ROBUSTNESS DATA FILE FOR THE DISCRETE ESTIMATOR DESIGNED 

AUTOPILOT 



FILE: ESXIMArO OmTA A I 

LOGO 0 



1C0.0 


0 


.0 


1.0 








1 




0 


0 








0.4 


0 


.4 


0.0 0 








0 0 0 




5 


3 i 1000 


0 






1717 170217 1 702 17 1 7 












1* 5 7957 58b- Gl 




-X 


. 6 3o9 32 6c — J 2 


5 • 53465 92E -03 


-3. 5* 1 5 2 3 cE - 2 


-5. 6350634E-3 


1. 4376905E- 03 




- 


o. 355CI95E-3 


3 - 12 1 50 3 80 E+ 0 


-2. 7434696E+0 


8 • *>960 341 c — 0 1 


0 

0 






0 

0 


0 


0 


0 


7. 37021215-03 




1 


.0001 347E-00 


1. 214571 3E -02 


-1. 31 10829E-4 


-2. 3707903E-5 


6.48032 77E-C6 




- 


*♦. j 148 jo 75- 5 


2 . J5 544 1 5E -02 
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